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ABSTRACT 



^* Received 

O . 

■ 

CO 

^"—•j Context. The open cluster NGC 1893, illuminating the Hn region IC 410, contains a moderately large population of O-type stars and is one of 
j — ' the youngest clusters observable in the optical range. It is suspected to harbour a large population of pre-main-sequence (PMS) stars. 
f"*) , Aims. We have probed the stellar population of NGC 1893 in an attempt to determine its size and extent. In particular, we look for signs of 
sequential star formation. 

Methods. We classify a large sample of cluster members with new intermediate resolution spectroscopy. We use Ha slitless spectroscopy of the 
field to search for emission line objects, identifying 18 emission-line PMS stars. We then combine existing optical photometry with 2MASS 



Oh 
i 



O ' JHKs photometry to detect stars with infrared excesses, finding close to 20 more PMS candidates. 



Results. While almost all stars earlier than B2 indicate standard reddening, all later cluster members show strong deviations from a standard 



C/5 reddening law, which we interpret in terms of infrared excess emission. Emission-line stars and IR-excess objects show the same spatial 
P*J distribution, concentrating around two localised areas, the immediate vicinity of the pennant nebulae Sim 129 and Sim 130 and the area close 
to the cluster core where the rim of the molecular cloud associated with IC 410 is illuminated by the nearby O-type stars. In and around the 
V~J ■ emission nebula Sim 130 we find three Herbig Be stars with spectral types in the Bl-4 range and several other fainter emission-line stars. 
' We obtain a complete census of B-type stars by combining Stromgren, Johnson and 2MASS photometry and find a deficit of intermediate 
^ mass stars compared to massive stars. We observe a relatively extended halo of massive stars surrounding the cluster without an accompanying 
population of intermediate-mass stars. 

Conclusions. Stars in NGC 1893 show strong indications of being extremely young. The pennant nebula Sim 130 is an area of active massive 
star formation, displaying very good evidence for triggering by the presence of nearby massive stars. The overall picture of star formation in 
NGC 1893 suggests a very complex process. 

Key words, open clusters and associations: individual: NGC 1893 - stars: pre-main sequence - stars: emission line, Be - stars: early-type 



1. Introduction 

High mass stars are known to form preferentially in star clus- 
ters, but the exact details of how they are born and whether 
their formation has an impact on the formation of less massive 
stars are still the arguments of open discussions (see, e.g., ref- 
erences in Crowther 2002). As massive stars disrupt the molec- 
ular clouds from which they are born, it is generally assumed 
that the formation of massive stars closes a particular star for- 
mation episode by eliminating the material from which further 
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* Partially based on observations obtained at the Nordic Optical 
Telescope and the Isaac Newton Telescope (La Palma, Spain) and 
Observatoire de Haute Provence (CNRS, France). 



stars may form (e.g.. Franco et al.ll 19941) . In this view, the for- 
mation of massive stars in an environment must take place over 
a short timescale. 

Numerous examples, however, seem to support the idea 
that the presence of massive stars triggers the formation of 
new stars in the areas immediately adjacent to their loca- 
tion (e.g., Walbornl 2002h . Such scenario would explain the 
form ation of OB associations, extending over dozens of par- 
sees dElmegreen & Ladalll977l) . Sequential star formation has 
been observed over both rather small sp atial scales (e.g., 
Deharveng et al . 2003; IZavagno et al . 2006) and large, massiv e 



stellar complexes, such as 30 Dor (Wa lborn & Blades 1997 ) 
but in most cases doubts arise about the role of the first gen- 
eration of stars: does their impact on the surrounding medium 
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actually trigger new star-formation episodes or simply blows 
away the clouds surrounding regions where star formation was 
already happening anyway? 

Investigations aimed at studying these q uestions can take a 
statistical approach (as in Massev et aljl995 ) or concentrate on 
the detailed study of one particular open cluster where star for- 
mati on is known to occur ( e.g., the investigation of NGC 661 1 
by Hillenbrand et al.lll993f) . Unfortunately, there are not many 
open clusters with active star formation and a large population 
of OB stars easily accessible for these studies. 

One such cluster is NGC 1893, a rather massive cluster, 
with five catalogued O-type stars, which seems to have very re- 
cently emerged from its parental molecular cloud. The ionising 
flux of the O-type stars has generated the H n region IC 410 on 
the edge of the molecular cloud (see Fig . [TJ . Though NGC 1893 
is rather more distant than some other areas where star forma- 
tion can be studied, its very young age, moderately rich O star 
population and relatively low interstellar reddening make it a 
very interesting target for optical/infrared studies. Moreover, its 
large Galactocentric distance and projection onto a molecular 
cloud mean that there is very little contamination by a back- 
ground population. 

As with all very young open clusters, the age of NGC 1893 
is uncertain. Several authors have estim ated it around ~ 4 Myr 
(Tapi aet al.lll991t IVallenari et al. Il999h . but the presence of 
luminosity class V early O-type stars would indicate an age 
< 3 Myr. Moreover, some emission-li ne B-type stars found 



in thi s area are li kely to b e PMS stars ( Marco & Neg ueruela 



20021 henceforth IPaper IS) , and hence very young (as con- 
traction time s for B-type stars are < 1 Myr). Even worse, 
Massev et alj (1 1995b classify some early-B stars in NGC 1893 
as luminosity class III, implying ages ~ 10 Myr. However, 
Massev et al] (119951) indicate that their spectral classifications 
are relatively rough, being directed to obtaining average spec- 
troscopic distances rather than discussing actual evolutionary 
stages. 

In Marco et al. henceforth IPaper il) . we used 

ubvy Hy6 CCD photometry of ~ 40 very likely main-sequence 
(MS) members to derive E (b - y) = 0.33 + 0.03 and V Q - M v 



13.9+0.2 for NGC 1893. In lPaper III we identified several PMS 
candidates, based on their spectral type and observed colours, 
three of which were shown to be emission-line PMS stars. 

In this paper, we investigate the possible age spread in 
NGC 1893 and take a fresh look at the star formation process 
in this area by considering a rather larger field. The paper is 
structured as follows: in Section [2] we present the new obser- 
vations used in this study, which we discuss in Section [3] to- 
gether with existing optical photometry. In Section @] we use 
2MASS JHKs photometry to find stars with infrared excesses, 
which we identify as PMS candidates. We show that candidates 
concentrate around only two locations, notably in the vicinity 
of the bright emission nebula Sim 130. In Section|5] we present 
a spectroscopic study of stars in the area of this nebula, identi- 
fying several emission-line stars. Finally, in Section|6] we dis- 
cuss the interpretation of the results in terms of evidence for 
triggered star formation. 



Table 1. Log of new spectroscopic observations. The top panel 
shows low resolution spectroscopic observations. The bottom 
panel displays intermediate resolution observations. February 
2001 observations are from the INT. October 2001 observations 
were taken with the OHP 1 .93-m. December 2001 observations 
were taken with the NOT. 



Star 


Date of 


Dispersion 


A Range 




observation 






S1R2N38 


2001 Dec 5 


1.5 A/pixel 


3830-6830 A 


S2R2N43 


2001 Dec 5 


1 .5 A/pixel 


3830-6830 A 


S1R2N26 


2001 Dec 6 


3. OA/pixel 


3100-9100 A 


S1R2N26 


2001 Dec 7 


2.3A/pixel 


3100-6675 A 


S5003 


2001 Oct 25 


1.8 A/pixel 


3800-6900 A 

JUVU V/V./V 1 V 


S5003 


2001 Dec 6 


3. OA/pixel 


3100-9100 A 


E09 


2001 Dec 6 


3. OA/pixel 


3100-9100 A 


E10 


2001 Dec 7 


3. OA/pixel 


3100-9100 A 


E17 


2001 Dec 7 


3. OA/pixel 


3100-9100 A 


^ 1K2IN 3 J 


2UU1 reb o 


0.4A/pixel 


oncn ennn a 

jyju-iuuu a 


C 1D1M1 A 

MK2JN 14 


2UU1 (JCt 24 


U.9A/pixel 


11 A C CCOC A 

j /4j— jd Ij A 


J> 1K2JN4U 


2UU1 Oct 2 j 


A C\ A In', .r ~1 

U.9A/pixel 


J /4J—JD/J A 




9001 rw 9 A 


u.7/\/pixei 


J /HJ — J J I J r\ 


S1R2N55 


2001 Oct 22 


0.5A/pixel 


6250-7140 A 


S1R2N55 


2001 Oct 23 


0.9A/pixel 


3745-5575 A 


S1R2N56 


2001 Oct 22 


0.5 A/pixel 


6250-7140 A 


S1R2N56 


2001 Oct 23 


0.9A/pixel 


3745-5575 A 


S1R3N35 


2001 Oct 23 


0.9A/pixel 


3745-5575 A 


S1R3N48 


2001 Oct 24 


0.9A/pixel 


3745-5575 A 


S2R3N35 


2001 Oct 24 


0.9A/pixel 


3745-5575 A 


S2R4N3 


2001 Oct 24 


0.9A/pixel 


3745-5575 A 


S3R1N5 


2001 Feb 9 


0.4A/pixel 


3950-5000 A 


S3R1N16 


2001 Feb 9 


0.4A/pixel 


3950-5000 A 


S3R2N15 


2001 Feb 9 


0.4A/pixel 


3950-5000 A 


S4R2N17 


2001 Feb 9 


0.4A/pixel 


3950-5000 A 


Hoag 7 


2001 Oct 23 


0.9A/pixel 


3745-5575 A 


HD 243035 


2001 Oct 24 


0.9A/pixel 


3745-5575 A 


HD 243070 


2001 Oct 24 


0.9A/pixel 


3745-5575 A 



2. Observations and data 

2.1. New data 

We obtained imaging and slitless spectroscopy of NGC 1893 
using the Andalucia Faint Object Spectrograph and Camera 
(ALFOSC) on the 2.6-m Nordic Optical Telescope (NOT) in 
La Palma, Spain, on the nights of December 5th-7th, 2001. 
The instrument was equipped with a thinned 2048 x 2048 
pixel Loral/Lesser CCD, covering a field of view of 6.'4 x 6.'4. 
Standard Bessell UBVRI filters were mounted on the filter 
wheel, while a narrow -band Ha filter (filter #21, centred on 
A = 6564A and with FWHM= 33A) was mounted on the FASU 
wheel. 

For the slitless spectroscopy, we made use of the Bessell 
R filter and grism #4. In total, 5 slightly overlapping images 
were taken, with 900-s exposure times. The weather was rel- 
atively poor, with some thin cloud veiling present on some of 
our exposures. The area covered by these observations is shown 
in Fig. Q] Spectroscopy of several emission-line stars in the 
field was taken on the same nights using the same instrument. 
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Fig. 1. The approximate boundaries of the five frames observed with slitless spectroscopy are indicated on the NGC 1893 field 
from a digitised DSS2-red plate. Note the position of the emission nebulae Sim 129 and Sim 130. The large white patch almost 
devoid of stars to the SW is the molecular cloud associated with IC 410. The five O-type stars in NGC 1893 are marked by circles 
and identified by name. The 8 early B-type stars in the periphery of NGC 1893 for which we present classification spectra are 
identified by squares. The O-type stars HD 242908 and HD 242926 are surrounded by bright nebulosity (darker patches). 



Grisms #3, #4 and #7 were used. A list of the objects observed 
is given in Table 1 . 

We obtained intermediate-resolution spectra of stars in the 
region of the bright nebula Sim 130 and surrounding area dur- 
ing 22nd-25th October 2001, using the 1.93-m telescope at 
the Observatoire de Haute Provence, France. The telescope 
was equipped with the long-slit spectrograph Carelec and the 
EEV CCD. On the night of 22nd October, we used the 1200 
ln/mm grating in first order, which gives a nominal disper- 
sion of * 0.45 A/pixel over the range 6245 -7 145 A. On the 
nights of 23rd & 24th October, the 600 ln/mm grating was used, 
giving a nominal dispersion of « 0.9A/pixel over the range 
3745-5575A. Finally, on 25th October, the 300 ln/mm grat- 
ing was used, giving nominal dispersion of * 1.8 A/pixel over 
the 3600 -6900A range. 



Finally, intermediate-resolution spectroscopy of several 
bright stars in the field of NGC 1893 was obtained during a 
run at the 2.5-m Isaac Newton Telescope (INT) in La Palma 
(Spain), in Febr uary 200 1 . Details on the configurations used 
can be found in Paper I j, while a list of all the observations 
presented here is given in Table 1 . 



All the data have bee n reduced using the Starlink 



software packages ccd pack ([Draper et al. 2000) and figaro 



(Shortridge et al. I l997l) and analysed using figaro and dipso 
(Ho warth et al.ll998l) . Sky subtraction was carried out by using 
the POLYSKY procedure, which fits a low-degree polynomial 
to points in two regions on each side of the spectrum. The ex- 
tent of these regions and their distance to the spectrum were 
selected in order to reduce the contamination due to nebular 
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emission. Bright sky lines coming from diffuse nebular emis- 
sion are visible in almost all of our spectra. 

2.2. Existing data 

For the analysis, we combine the new observations with exist- 
ing photometric datasets. On one side, we use JHK S ph otom- 
etry from the 2MASS catalogue dSkrutskie et al.ll2006l) . The 
completeness limit of this catalogue is set at Ks, - 14.3, which 
- as we shall see - roughly corresponds to the magnitude of an 
early A-type star in NGC 1893. In addition , we have two opti- 
cal photometric datasets. UBV data from lMassey et al.1 (11995b 
cover a wide area around the cluster. Its magnitude limit is 
close to V ~ 18, but errors become a ppreciable for V > 16. 
Stromgren photometry from Paper il only reaches V ~ 16 
and covers a more restricted central area. However, Stromgren 
photometry allows the determination of approximate spectral 
types un der the single assumption of standard reddening. From 
Paper fl, we know that V ~ 16 roughly (depending on the red- 
dening) corresponds to the magnitude of an early A-type star in 
the cluster. Therefore, pretty much by chance, all three datasets 
have comparable limits. 



3. Results 

3.1. Emission-line stars 

In this paper, we address the search of emission-line PMS stars 
following a new approach: deep slitless spectroscopy of the 
whole field. This technique, based on the use of a low disper- 
sion grism coupled with a broad-band filter resulting in an "ob- 
jective prism-like " spectrogram of all the o bjects in the field, 
has been used by iBernabei & Polcarol d200ll) for the search of 
emission line stars in open clusters. By combining a Johnson R 
filter and a low-resolution grism, we obtain bandpass imaging 
spectroscopy centred on Ho-. 



Table 2. Known emission-line stars in NGC 1893. E01 and E02 
w ere alrea dy listed in the literature. E03-E07 were described 
in Paper III E08-E18 are found or confirmed here. The spec- 



The obvious advantage of this technique with respect to 
narrow-band photometry (as used in lPaper II) is that it can reach 
very faint stars. The detection limit is difficult to define. In un- 
crowded regions, it is reached when the spectra are too faint 
to be seen against the background, but in crowded regions the 
overlap of adjacent spectra becomes im portant. Comparison 
to the photometry of iMassey et al.1 (11995b shows that we have 
been able to detect emission lines in stars fainter than their limit 
at V ^ 18, but we cannot claim completeness. 

We detect t he 7 pre viously known emission-line objects 
(the 5 listed in iPaper III and the two catalogued Ho- emitters 
close to Si m 130). W e fail to detect the candidate emission- 
line S5003 (IPaper HP . as its spectrum, given the orientation of 
our dispersion direction, is completely superimposed by that 
of its brighter neighbour S3R1N9. However, we have obtained 
a long-slit spectrum of this object and can confirm it as an 
emission-line star (see Table 13. U . In addition, we detect 10 
new stars with Her emission (see Table 13. U . and name them 
E09 - El 8. Only two of the new emission-line objects are bright 
enough to have been observed by previous authors, namely 
S1R2N55 and S1R2N26. We confirmed the emission-line na- 
ture of some of these objects through long-slit spectroscopy 



tral types of E09, Ell and E17 cannot be determined from our 
spectra. We did not take long-slit spectra of objects whose spec- 
tral type is marked as Note that the Ks magnitudes (all 
from 2MASS) of some objects are affected by blending. 



Name 


RA 


Dec 


Spectral 












Tvnp 




E01 


=S1R2N35 


05 23 09.2 


+33 30 02 


B1.5 Ve 


10.01 


E02 


=S1R2N38 


05 23 04.3 


+33 28 46 


B4Ve 


10.25 


E03 


=S3R1N3 


05 22 43.0 


+33 25 05 


B0.5 IVe 


10.9 


E04 


=S3R1N4 


05 22 46.1 


+33 24 57 


B1.5 Ve 


12.3? 


E05 


=S2R1N26 


05 22 48.2 


+33 25 00 


~G0Ve 


11.9? 


E06 


=S2R1N16 


05 22 51.1 


+33 25 47 


~F7Ve 


11.8 


E07 


=S1R2N23 


05 22 52.1 


+33 30 00 


~F6Ve 


11.3 


E08 


= S5003 


05 22 40.8 


+33 24 39 


Ke 


13.5 


E09 




05 22 43.8 


+33 25 26 


? 


9.4 


E10 




05 22 49.6 


+33 30 00 


? 


>15 


Ell 


= S1R2N26 


05 22 57.9 


+33 30 42 


-A3 Ve 


12.6 


E12 




05 23 00.0 


+ 33 30 41 


? 


13.0 


E13 




05 23 02.8 


+ 33 29 40 




13.9 


E14 




05 23 04.4 


+ 33 29 48 




13.2 


E15 




05 23 06.3 


+33 31 02 




13.2 


E16 


= S1R2N55N 


05 23 08.3 


+33 28 38 


B1.5 Ve 


> 10.3 


E17 




05 23 08.9 


+ 33 28 32 


? 


11.5 


E18 




05 23 09.9 


+ 33 29 09 




13.2 



(Table 1). Most of them are too faint for spectral classifica- 
tion. S1R2N26 is an A-type star with strong Ha emission. E12 
does not show any photospheric features, but has all Balmer 
lines in emission and also shows strong Ca n K lines and Ca n 
8498, 8542,8662 A triplet emission. All the known emission- 
line stars in NGC 1893 are listed in Table O 



3.2. O-type stars 

We have obtained accurate classifications for four O-type stars 
in the field of NGC 1893. The fifth O-type star, to the South of 
the cluster (HD 242926), was not observed here, but it has been 
observed as part of another pr ogramme and it s spectral type is 
07 V, in total agreement with IWalbornl dl973 ). This star is re - 
ported to show strong radial velocity changes by|j ones ( 1972). 

The spectra of the four O-type stars observed are pre- 
sented in Fig. [2] The spectral types have bee n derived using the 
quantitative methodology of Mat hvsl dl988l) . based on Conti's 
scheme. For S4R2N17 (HD 242908), the spectral type criteria 
fall on the border between 04 a nd Q5 and we will adopt the 
04V((f)) classification given by lWalbornldl973h . In S3R2N15 
(LS V +34° 15), the He i lin es are stronge r and the quantitative 
criteria indicate 05.5 V((f)). |Jonesldl972l) reports two measure- 
ments of the radial velocity of this star, differing by more than 
100 km s _1 , strongly suggesting that there are two O-type stars 
in this system. In S3R1N16 (BD +33° 1025 A), the condition 
He i /14471A =a Hen A4541A implies by definition 07, while 
the strength of He n /I4686A and very weak N m emission in- 
dicate a MS classification. Finally, in S3R1N5 (HD 242935), 
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Fig. 2. Classification spectra of the 4 O-type stars near the core 
of NGC 1893. The earliest spectral type is that of S4R2N17 
(HD 242908), 04 V((f)). The other three objects have spec- 
tral types 05.5V((f)) for S3R2N15, 07 V for S3R1N16 and 
07.5V((f))forS3RlN5. 



He i 14471 A is slightly stronger than Hen ,14541 A and the 
quantitative criteria indicate 07.5 V((f)), though Hen ,14686 A 
is rather weak and clos e to the limit for luminosity class III 
given bv lMathvsl dl988l) . 



3.3. B-type stars 

We have also derived new classifications for bright B-type 
stars to the East of NGC 1893, in order to check if they 
are connected to the cluster in spite of their relatively large 
angular distance to the cluster core (see their distribution in 
Fig. [TJ. These spectra are displayed in Fig. [3] At the resolution 
of our spectra, the traditional criterion for spectral classifica- 
tion around B 0, namely the ratio between Si iv /14089A and 
Sim 14552A dWalborn & Fitzpatrick|[l990l) is difficult to ap- 
ply, since Si iv /14089A is blended into the blue wing of H6. 
We have thus resorted to additional criteria. 

S1R3N48 is the earliest object in the sample, and the only 
one for which we give a B0V classification, based on the 
conditions Hen 14686A^Cm ,i4650A and Hen ,14686A> 
He i /14713A. We have classified as B0. 2 V those stars in which 
Hen /14200A is barely visible and Hen 14686A is compara- 
ble in strength to He i /14713A. We have taken as B0. 3 V those 
stars in which He n /14200A is not visible and He n 14686A is 
clearly weaker than He i ,147 1 3 A, and assigned B0. 5 V to those 
objects in which Hen /14686A is visible, but very weak. All 
the stars analysed fall in the B0 -B0. 5 range (see Table 3), ex- 
cept for S2R3N35. For this object, we derive a spectral type 
B2.5 V, In this spectral range, the luminosity class is mainly in- 
dicated by the strength of the Si in and Si iv lines compared to 
those of He i, while C m and O n lines can be used as subsidiary 
indicators (see Walborn & Fitzpatricklll990l for a detailed dis- 
cussion). None of the spectra analysed suggests a luminosity 
class different from V. 



Table 3. O-type and early B-type stars around NGC 1893 
for which we derive accurate spectral types, tog ether with the 
derive d reddening. Photometric data are from iMassey et al 
(1995). For HD 242935 there are no reliable photometric mea- 
surements due to heavy blending. 



Star 


Name 


V 


Spectral 


E(B - V) 


Number 






Type 




S4R2N17 


HD 242908 


9.03 


04V((f)) 


0.57 


S3R2N15 


LS V+34°15 


10.17 


05.5 V((f)) 


0.79 


S3R1N16 


BD +33°1025A 


10.38 


07 V 


0.59 


S3R1N5 


HD 242935 




07.5 V((f)) 




S1R2N14 


LS V +33°23 


11.22 


B0.2V 


0.45 


S1R3N35 


LS V +33°26 


11.17 


B0.2V 


0.54 


S1R3N48 


HD 243018 


10.94 


B0V 


0.42 


S2R3N35 




12.43 


B2.5V 


0.47 


S2R4N3 


LS V +33°24 


11.02 


B0.5V 


0.66 




HD 243035 


10.90 


B0.3V 


0.50 




HD 243070 


10.83 


B0.2V 


0.51 


Hoag 7 


LS V +33°27 


10.74 


B0.3V 


0.49 



For all the B -type st ars with accurate spectral types ei- 
ther here or in IPaper it we ca lculate the d i stance modu- 
lus, using UB V magnitudes f rom iMassey et al. ( 1995b . intrin- 
sic colours of iWegnerld 1994b and absolute magnitudes from 
iHumphrevs & McElrovl dl984b . under the assumption of stan- 
dard reddening (justified in the next section). The average value 
is 13.4 ± 0.3, in g ood agreement with the value found by 



Massev et ail d 1995b . 



4. 2MASS data 

4.1. Interstellar Reddening 

The reddening along the face of NGC 1893 is known to be 
variable and this may have a bearing on the derivation of clus- 
ter parameters, as most methods are very sensitive to the in- 
terstellar reddening law assumed and treatment of the redden- 
ing. Deviations from the standard value are frequent in the 
optical, though the reddening law in the infrared has been 
prove d to show very lit t le va riability along different lines of 



sight (Indebe touw et al.l 120051) . Because of this, we use the 



2MASS JHKs photometry to check if the extinction law to- 
wards NGC 1893 is standard. 

For all B-type member stars with accurate Stromgren pho- 
tometry in IPaper J, we calculate individual values of E(B - V) 
by using the simple relation E(B — V) = \AE(b - y). From 
the 2MASS magnitudes, we calculate E(J - Ks) assuming 
the intrinsic colours from iDucati et al.l (2001) and the spec- 
tral types derived fr om the Stromgren photo metry. For a stan- 
dard reddening law ( Rieke & Lebofskyl 1985b . we should have 
E(J - K s ) - 0.5E(B - V)Q. We find that many stars have 
E(J - Ks) ~ 0.5E(B - V), but a significant fraction show rather 
larger E(J - Ks) than expected from their E(B - V) . 



1 The extinction law o f lRieke & Lebofskvl |l985) uses Johnson's K 
instead of i^ s , but the difference is negligible. Using the values calcu- 
lated by Hanson (2003), we would have E(J - A* s ) = 0.48£(B - V) 
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Fig. 3. Classification spectra of B-type stars in the area of NGC 
1893. From top to bottom, stars are displayed from earlier to 
later spectral type. Photospheric features used for spectral clas- 
sification in this range are indicated. 



This deviation from the standard relationship is very un- 
likely due to a non-standard extinction law, because the amount 
of excess E{J - Ks) is not strongly correlated with position 
within cluster. On the other hand, Fig. H] shows the dependence 
of E( J - Ks) with spectral type. For all O and B-type members, 
we plot E(J - Ks) against the T e s corresponding to the spec 



tral ty pe derived according t o the calibrations o f Martins et al. 
d2005l) for O-type stars and lHumphrevs & McElrovl dl984l) for 



B-type stars. Except for three stars, all stars earlier than B3 
have about the same E(J - Ks), compatible or very slightly 
above the value expected for a standard reddening lawQ All 
stars with spectral type B3 or later (except for two) have larger 
excesses, with a clear tendency to have larger excesses as we 
move to later spectral types. This dependence of the infrared 
excess on spectral type, while it shows no dependence on loca- 
tion, is clearly suggesting that the excesses are intrinsic to the 
stars and not related to the extinction law. 

To investigate this further, we use the chorizos code 
( Maf z-Apellania 120041) to estimate the extinction law. This 




Fig. 4. Plot of the infrared excess E(J - Ks) against spectral 
type (represented by the corresponding T e g) for O and B-type 
MS stars in the cluster. The straight line E(J - K s) = 0.26 cor- 
respo nds to the cluster average E(B — V) = 0.53 dMassev et al 



2 Note that S3R2N15 has a high E(J - K s ), but it also has E(B - V) 
much above other stars in Table 3. 



19951) if a standard reddening is assumed. The error bars in 
E(J - Ks) represent only the photometric errors. The uncer- 
tainty in the intrinsic colours has not been included, as it is 
difficult to quantify and unlikely to depend on spectral type. 



program tries to reproduce an observed energy spectral dis- 
tribution by fitting extinction laws from ICardelli et al.l 0989) 
to the spectral distribution of a stellar model . As input, we 
used the UBV photometry from iMassev et al.l d 19951) and the 
JHKs photometry from 2MASS, together with Kurucz models 
of main sequence stars with the r e jf corresponding to our stars 
and logg = 4.0. We run the program for all the stars with spec- 
tral types derived from spectra or from Stromgren photometry. 
Out of 20 stars earlier than B3, 15 are best fit by reddening 
laws having 2.8 < R < 3.4, while 5 require R > 3.5. Out of 28 
stars later than B3, 25 require R > 3.5 and 3 did not converge 
to a solution. Again, we find a strong dependence of the red- 
dening law on the spectral type. Taking the average for all the 
stars with spectral type earlier than B3, we find R = 3.3 ± 0.2. 
Leaving out the 5 stars with R > 3.5, we find an average 
R = 3.16 + 0.12. 

Our interpretation of this result is that the interstellar red- 
dening law to NGC 1893 is standard, but a substantial num- 
ber of objects show important individual (J - Ks) excesses. 
These excesses are present in all stars later than B3 and in a 
few early stars (these early-type stars with anomalous values 
of R may perhaps have later-type companions with individ- 
ual excesses). This interpretation is further supported by the 
Ks/(J - Ks) diagram for cluster members (Fig. |5). All stars 
of mid and late B spectral type deviate strongly from the al- 
most vertical main sequence traced by early members, display- 
ing much larger (J - Ks). This separation is in clear contrast 
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Fig. 5. Plot of K s against (7 - # s ) for all the 2MASS stars ful- 
filling the two conditions set in Section l4~2l Fill ed circles r ep- 
resent O and B-type MS members from lMarco et al.l (1200 II) or 
with spectra in this paper. Open circles are known emission- 
line PMS stars. Stars represent the rest of the IR excess can- 
didates, three of which fall within the main sequence traced 
by known members (one has a known member superimposed). 
Note also the object with (J-Ks) =s 3, our emission-line source 
E09. 

with the fact th at stars in the B0- B8 range have almost iden- 
tical (J - K)q dDucati et al.l 120011) and again confirms that all 
stars later than B3 show some E(J - Ks) excess, with stronger 
excesses corresponding to later spectral types. 

Mid and late B stars fit the ZAMS rather well in differ- 



ent optical observa tional HR diagrams (e.g., Tapi a et al.111991 



Masse v et al. 1995; Marc o et al ] |200ll) . though the y do not oc- 



cupy standard positions in the [ci]/[mi] diagram dTapia et al 



1 99 it iMarco et al.ll200 II) . How do we then explain their E(J 



Ks) excesses? In Fig. [6] we plot an infrared HR diagram for the 
B-type members. For each star, we assume that the interstellar 
reddening E(J- Ks) ls = Q.5E(B-V) and the corresponding ex- 
tinction is Ak s = 0.67£ , (7-/Ts)i s , according to the standard red- 
dening law. We then shift their = K$ - Ak s by DM= 13.5 
(see Section 16.1b and plot them in the colo ur-magnitude dia- 
gram. We also plot the PMS isochrones from lSiess et al.l d2000) 
for 1 and 2 Myr, t ogether with the youngest isochrone from 



Girardi et al.1 (12002), corresponding to 4 Myr, as the position of 



most B stars on this isochrone will not deviate in any measur- 
able way from the ZAMS. 

We observe that the smooth way in which the E(J - Ks) 
excesses increase with decreasing mass does not fit at all the 
shape of the PMS isochrones. As a matter of fact, the posi- 
tions of the stars in the HR diagram cannot mean that they are 
moving towards the main sequence along PMS tracks, because 
the infrared excess, as it has been defined, is a measure of the 
discrepancy between optical and infrared colours. PMS stars 



t-10 yr 




(J-K,) 



Fig. 6. Observational HR diagram for B-type members of 
NGC 1893. The loc ation of the ZAMS i s represented by the 4- 
Myr isochrone from Girardi et al. (2002). The PMS isochrones 
for 1 and 2 Myr from lSiess et al. (2000) are also shown. 



should have both optical and infrared colours appropriate for 
the position in the theoretical HR diagram. Indeed, if the devi- 
ation of B3 - 4 stars from the ZAMS had to be attributed to their 
being on PMS tracks, the age of the cluster should be only a few 
10 5 yr. In view of this and the fact that most B-type stars fit the 
ZAMS well in the optical, we are led to interpret the E(J - Ks) 
excesses as due to the presence of material left over from the 
star formation process, most likely in the form of a remnant 
disk. If so, the PMS isochrones would suggest an age for the 
cluster of ~ 2 Myr, with most stars later than B2 (i.e., stars with 
M* < 8Mq) showing evidence for some remnant of a disk. This 
remnant cannot be very important, because stars fit the ZAMS 
well in the optical (both Johnson and Stromgren) HR diagrams 
and are not detected as emission line objects. High SNR spec- 
tra of the Ha line or in the H or K bands may be able to reveal 
some spectroscopic signature of such remnants and so confirm 
this hypothesis. 

4.2. OB stars and PMS selection 

The OB stars observed extend over a wide region to the East 
of the molecular cloud associated with IC 410. The core of 
NGC 1893 (defined as the region with the highest concentration 
of MS members) lies immediately to the West of the molecular 
cloud. There is an abrupt decrease in the number of optically 
visible stars as we move East from the O-type stars S3R1N16, 
HD 242935 and S3R2N15 (see Fig.Q). This suggests that an 
opaque part of the molecular cloud blocks our view (see also 



Leisawitz et al. 19891) . As we lack three dimensional informa- 



tion, we cannot know if the cluster is spread over the wall of 
the dark cloud, but the fact that there is no bright luminosity 
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in this area suggests that the dark cloud is partially between us 
and the cluster. As large areas of dark nebulosity block some 
sight-lines, we consider the possibility that there may be other 
OB stars in the field, obscured by gas and dust. 

The /, H and Ks magnitudes from the 2MASS catalogue 
can be used to look for reddened early-type stars, under the 
assumption of a standard reddening la w. Taking Mk — -1 .6 as 
the intrinsic magnitude fo r a B2 V star (Humphreys & McElroy 
1984 bucati et alj|200lh and DM = 13.9 to NGC 1893 (con- 



sidered an upper limit), all stars earlier than B2 located at the 
distance of the cluster an d reddened according to the law of 



Rieke & Le bofsky ( 1985) will fulfil the following condition: 



K s - 1.78(7/ -tfs) < 12.5 



(1) 



Obviously many other stars will also fulfil this condition. 
However, if we define the reddening free parameter Q — (J - 
H)-UO(H-Ks), OB stars will have Q =* 0.0. Combination of 
Condition [T] with Q 0.0 has been sho wn to be very efficient 
at ide ntifying reddened OB stars (e.g., IComeron & Pasq uali 
20051) . However, not only OB stars fulfil both conditions. 



Foreground A-type stars and background red giants and super- 
giants will also fulfil the conditions. However, in our case, con- 
tamination by background red stars is unlikely, as the cluster 
lies at a large Galactocentric distance in the Anticentre direc- 
tion and is projected on to a dark cloud. 

We take all 2MASS objects within a 15' radius of the centre 
of the cluster with magnitudes tagged as good and an error in 
the Ks magnitude 6Ks < 0.05 mag (larger errors would imply 
completely unreliable Q parameters) and select all stars ful- 
filling Condition [TJ and having Q < 0.1. We discard foreground 
A-type stars using two criteria: (1) (B-K$) w indicates unred- 



dened stars (most stars have B magnitudes from lMassev et al.l 
d 19951) : for the rest we use USNO B1.0 magnitudes), and (2) 
(J - Ks) lower than the average for the known OB members 
indicates that the stars are foreground to the cluster (only ob- 
jects substantially more reddened than this average would have 
escaped optical surveys). 

All known members earlier than B2 are selected by these 
two criteria. Moreover, because of their infrared excesses, sev- 
eral B-type members later than B2 are also selected. As a 
matter of fact, a substantial fraction of the emission-line stars 
listed in Table 3 (the brightest ones), including some of F type, 
have been selected. This is not surprising as these stars have 
strong infrared excesses and then : (a) are much brighter in K$ 
than normal stars of the same spectral type (so much brighter 
that they pass our magnitude cut) and (b) the assumption of a 
standard interstellar law in the calculation of the Q parameter 
means that their (H - Ks) excesses are "over-corrected", re- 
sulting in negative values of Q, values that no normal star can 
have. 

In view of this, we take stars fulfilling Condition [TJ hav- 
ing Q < -0.05 and large values of (J - Ks) as infrared ex- 
cess objects and therefore PMS star candidates. This procedure 
does not select all the objects with infrared excess in the field, 
but only those relatively bright and with strong (H - Ks) ex- 
cesses. A search for a complete sample of infrared excess ob- 
jects in this field would need deeper photometry than provided 
by 2MASS and is beyond the scope of this paper. 



We find more than thirty stars fulfilling those criteria, 
among them, 11 of the known emission-line stars. Figure [5] 
shows the Ks/{J - Ks) diagra m for all MS member s (filled cir- 
cles) in NGC 1893 selected in lMarco et ail d200ll) and all the 
infrared excess stars. Known emission-line PMS stars (open 
circles) are located to the right of the main sequence, as ex- 
pected. We see that all the newly selected candidates occupy 
the same locations as the emission-line stars, but at fainter mag- 
nitudes, except for three objects falling close to the location of 
B-type MS members. Of t hese, on e is a photometric member 
outside the area covered in lPaper 1 according to its UBV mag- 
nitudes. A second one is very far away from the cluster and 
could only be a member if it is much more reddened than any 
other member. Hence we reject it as a good photometric can- 
didate. T he third object is S3R1N13, which was identified in 



Paper III as a candidate PMS star without emission lines. It has 



a spectral type B5 III-IV and its high reddening strongly sug- 
gests it is a PMS star in NGC 1893. The only alternative would 
be a foreground star with very unusual colours, but its spectrum 
does not show significant anomalies. 

Fig.|7]compares the distribution of these candidate infrared 
excess objects to that of MS members. The distribution is cer- 
tainly not random, as they strongly concentrate around two 
small areas, the vicinity of the pennant nebulae Sim 129 and 
Sim 130 and the rim of the molecular cloud closest to the clus- 
ter core. This distribution confirms beyond doubt that the stars 
selected are a population associated with the cluster rather than 
red background stars. Moreover, the spatial distribution of these 
objects coincides exactly with that of emission-line stars, giv- 
ing full support to the interpretation that most of them are also 
PMS stars. We cannot rule out the possibility that a few of 
these objects (especially the few ones at large distances from 
the cluster) are background red stars, but certainly the majority 
of these objects are PMS members with large infrared excesses. 

Once we account for known members, foreground stars, 
emission-line stars and infrared-excess candidate PMS stars, 
we have exhausted the list of stars selected according to our 
original criteria. This means that there are no obscured OB stars 
in this region, at least within the magnitude limit of 2MASS. 
Note, however, the position of E09 in Fig. [5] If this object be- 
longs to the cluster, it should be a young massive stellar ob- 
ject. Unfortunately, this object is so faint in the optical that 
our spectrum is extremely noisy. The only obvious feature is 
a very strong Ha emission line. Another possible emission fea- 
ture may correspond to the O i 8446A line. 

5. The area around Sim 130 

Results in the previous sections clearly show that present-day 
star formation in NGC 1893 is strongly concentrated towards 
the pennant nebulae Sim 129 and Sim 130, with most of the 
PMS stars located around the latter. 

The "head" of Sim 130 cont ains a group of stars which 



roup 



were observed photometrically by lTapia et a l. ( 1991) as if they 
were a single object (their Star 35). They derive the colours of 
an early-type star with emission lines. The head of Sim 130 
has also been identified as the near-infrared counterpart of the 
IRAS source IRAS 05198+3325, considered a Young Stellar 
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Fig. 7. The spatial distribution of all known members of NGC 1893 is shown on this DSS2 image. Objects on the main sequence 
according to their optical colours are shown as squares (red for spectral types B2 and earlier, yellow for B3 and later). PMS 
candidates (emission line stars and infrared excess objects) are shown as blue circles. Main-sequence objects are spread over a 
much larger area than PMS candidates. 



Object (YSO) candidate (CPM16 in lCampbell e"tal]|l989h . In 
addition, in the immediate vicinity of Sim 130, we find the 
emission-line star S1R2N35, which based on low-resolution 
spectra, was cons idered to be a very good candidate to a 
Herbig Be star in lPaper j . Not very far away lies the catalogued 
emission-line B star NX Aur = S1R2N38. 

Figure [8] identifies the brightest stars immersed in 
the bright nebulosity of Sim 130: S 1R2N56 ([MJD95 1 
J052307.57+332837.9 in the catalogue of lMassey et aD 19951) . 
S1R2N55 ([MJD95] J052308.30+332837.5) and a fainter 
star not observed by ICuffevI dl973l) . which we will call 
S5004 ([MJD95] J052306.7 1+332840.3). There are many 
other fainter stars within the nebula (more clearly seen in R- 
band images), among them our emission-line objects E17 and 
El 8. Also, partially immersed in the nebulosity, we find the 
bright starS!R2N44. 



5.1. S1R2N44 

S1R2N44 lies very close to Sim 130, just to the SW. Its spec- 
trum, displayed in Fig. [10] shows nebular emission lines on top 
of a normal absorption B-type stellar spectrum. A lower reso- 
lution spectrum does not show any signs of intrinsic emission 
in Ha. However, a cut of the spectrum in Ho- clearly shows 
that the intensity of nebular Ha emission increases consider- 
ably around the source, strongly suggesting the association of 
S1R2N44 with the nebulosity. 

From the Full Width at Half Maximum (FWHM) of four 
He i lines, we estimate for S 1R2N44 an apparent rotational ve- 
locit y v sin i «g 320 km s -1 , following the procedure described 
bv lSteeleetaII(ll999h. We estim ate ts spectral type at B2. 5 V. 



Fitzsimmonsl ( 1993 ) gives (b-y) = 0.273, imply- 



For this star. 

ing E(b-y) - 0.39, slightly above the average for cluster mem- 
bers. If the reddening is standard, V = 13.3 implies My = —1.9, 
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Fig. 8. An Ha image of the head of the cometary nebula Sim 
130 (a portion of a 300-s exposure obtained on the night of 5th 
December 2001 with ALFOSC on the NOT, equipped with the 
narrow-band Ha filter #21). The stars discussed in Section [5] 
have been identified. Note the bow-shaped emission rim and 
the extended nebulosity surrounding the stars. 
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Fig. 9. Ha spectra of the two stars in the "head" of Sim 130. 
While the emission features on the spectrum of S 1R2N56 seem 
to be entirely due to the surrounding nebulosity, S1R2N55 is 
clearly an emission-line star (see the inset for details). This ob- 
ject is the counterpart to the IRAS source 05198+3325, identi- 
fied as a massive young stellar object. 



which is not surprising for the spectral type. This star must have 
reached the ZAMS, as its observed (J - Ks) = 0.25 indicates 
that this object has no infrared excess. 

5.2. S1R2N56 

This star forms the "head" of Sim 130. Its spectrum is strongly 
contaminated by nebular emission, but this can be easily dis- 
tinguished from photospheric features at our resolution. There 
does not seem to be emission intrinsic to the star (see the Ho- 
profile in Fig. [9] where also a weak Hei /16678A absorption 
line can be seen). 

The blue spectrum of this object is displayed in Fig. \10\ 
From the FWHMs of four He i lines, we estimate an apparent 
rotational velocity v sin i « 2 10k m s" 1 . We estimate it s spectral 
type at B 1 .5 V. For this object, iMassey et al. (1995) measure 
(B -V) - 0.43, implying a reddening E(B - V) = 0.66, well 
above the average for cluster members. If the reddening law is 
standard, then the measured V = 13.54 implies My = -2.0, 
which, though slightly too low for the spectral type, is not 
in strong disagreement. The 2MASS colour {J - Ks) = 0.76 
implies substantial reddening and this object is selected as 
an infrared excess candidate, but this could mainly be due to 
contamination of its photometry by the bright nebulosity. The 
available evidence suggests that S1R2N56 is a star settling on 
to the ZAMS. 

Observations by previous authors seem to indicate very 
large variability in V with a long-term dimming from V < 13, 
but since none of these authors mentions S1R2N55, it is pos- 
sible that both stars have been measured together (this is cer- 
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Fig. 10. Blue spectra of S1R2N56 (top) and S1R2N44. 
S1R2N56 is immersed in the pennant nebula Sim 130 and 
its spectrum shows strong nebular emission. Nebular emis- 
sion is also present, though rather weaker, in the spectrum of 
S1R2N44, which lies just outside Sim 130. 



tainly the case in iTapia etalJll99ll whose star #35 corresponds 
toSlR2N56+SlR2N55). 

5.3. S1R2N55 

S1R2N55 is clearly immersed in the nebulosity associated with 
Sim 130. As a matter of fact, our NOT images clearly show 
S1R2N55 to be a close double, but the two components were 
not resolved during the OHP run, when spectra were taken. 

The red spectrum is shown in Fig. [9] In addition to strong 
nebular lines, a very broad Ha emission line can be seen: 
S1R2N55 is a Be star. The Ha line peaks at 15 times the con- 
tinuum intensity and has an Equivalent Width (EW) of -54 ± 3 
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Fig. 11. Blue spectra S2R1N35 (top) and S1R2N55, two early 
Herbig Be stars associated with Sim 130. S1R2N55 is actually 
a blend of two very close B-type stars, and likely only one of 
them is a Herbig Be star. 



A, a fraction of which is attributable to the nebular component. 
Given its presence in the middle of bright nebulosity in an H n 
region with active star formation, S1R2N55 is almost certainly 
a PMS Herbig Be star. In the /?-band, the Northern component 
is slightly brighter than the Southern one, but in the Ha im- 
ages, it is much brighter. This clearly shows that the Northern 
co mponent of S1R2N55 is the emission-line object. 

Ishii et all d200lh identified S1R2N55 as the counterpart to 



the IRAS source IRAS 05198+3325 (YSO CPM16). In their K- 
band spectrum, they found strong Bry emission from the source 
and H2 emission of nebular origin. 

The blue spectrum of S1R2N55 is shown in Fig. QT| 
Emission is present in HB and Hy. The spectrum is likely dom- 
inated by the brighter Northern component, but must be a com- 
bination of the spectra of the two stars. From the apparent pres- 
ence of weak 11 4076A and Si iv 4089A, the brighter compo- 
nent must be around B 1 V. The fainter component cannot have 
a very different spectral type, around B2 V. This object is se- 
lected among our infrared excess candidates. At least the bright 
component is a Herbig Be star. 

5.4. S5004 

The faint star S5004 is located just on the bright rim of nebu- 
losity defining the head of the cometary nebula. In spite of this, 
we have been able to clean the spectrum of nebular emission. 
The spectrum shows a prominent G-band, comparable in inten- 
sity to Hy, indicating a spectral type close to GO. The weakness 
of the Sr 11 4077 A and ot her luminosity indi cators seems com- 
patible with a GO V star. iMassey et al. I (1 1995b give V = 15.79, 
(B - V) = 1 .27 for this object. This star is prominent in the K- 
band (Ks = 11 .55) and is selected as one of our infrared excess 
candidates. Indeed, comparison of the observed (J -Ks) = 1.4 



ber, then its intrinsic magnitude would be My ~ —1.7, clearly 
very bright for a normal main-sequence G-type star. S5004 may 
then be an interme diate mass star still on t he contraction track, 
similar to G l Ori E ( Herbig & Griffinl2006l) . A spectrum around 
the Li 1 6707 A line could test this hypothesis. Alternatively, the 
photometry could be in error because of nebular contamination 
and this could be a foreground star. 

5.5. E17 

This star, too faint in the optical to have been observed by 
Massev et all d 1995b . is very bright in the ZT-band (K s = 1 1 .49 
in 2MASS) and is selected as an infrared excess object. Its 
spectrum is completely featureless, except for a prominent Ha 
emission line EW= -40 ± 2A. The emission line and infrared 
excess identify this object as a PMS star. 

5.6. S1R2N35 = E01 



with the colours expected for this spectral type (Ducati et al 
2001 ) implies E(J -Ks) = 1.13. This value is much larger than 
would correspond to its E(B - V), indicating a large infrared 
excess, typical of a PMS star. If this is object is a cluster mem- 



Classified by Kohoute k & Wehmeyer (1999) as an emission 
line star, this object in the immediate vicinity of Sim 130 was 
proposed as a Herbig Be candidate in lPaper j. based on a low- 
resolution spectrum. A higher resolution spectrum is shown in 
Fig. QT| Classification is complicated by the strong emission 
lines, but, based on the presence of some weak Oi lines and 
possible strength of Si m lines, while the Si iv lines are not vis- 
ibl e, we adopt a B 1 V, th ough it could be slightly later. 

IMassey et al] dl995h give V = 12.33, E(B - V) = 0.52, 
indicating a very large colour excess E(B - V) - 0.8. Again, 
this object is selected as an infrared excess candidate. From the 
intrinsic colour of a Bl V star, the implied excess is E(J-K)s = 
0.57. This object is hence a Herbig Be star. 

5.7. S1R2N38 = E02 

Already known as an emission-line and variable star (NX Aur), 
this object lies in the immediate vicinity of Sim 130, display- 
ing Ha and HB strongly in emission, EW(Ha)= -55 + 2A and 
EW(Hy8)= -3.7 + 0.3A. The upper Balmer lines display weak 
blue-shifted emission components. Prominent emission lines 
of Fen are also present. Though the object is clearly a B-type 
star, an exact spectral type is difficult to derive from our low- 
resolution spectrum. Based on the strength of the Mgi 4481 A 
lin e, we estimate i t to be B4 V. 

Massey et al.1 d 19951) give V = 14.41, E(B - V) = 0.57, 
indicating a E(B - V) « 0.8. 2MASS gives K s = 10.25 and 
(J - Ks) = 1.76, implying a huge infrared excess E(J - K$) ~ 
1 .9. This object is hence also a Herbig Be star. 

6. Discussion 

6.1. Distance, reddening and extent 

Our analysis shows that the optical/near-IR spectral energy dis- 
tributions of almost all stars earlier than B3 are best fit when a 
standard R = 3.1 reddening law is assumed. A s t andar d inter- 
stellar law has also been found by Yadav & Sagar ( 2001), using 
less sophisticated techniques. Tapia et alj dl991 ) found a value 
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of R — 2.8 ±0.1, slightly lower than our value. On the other 
hand, all stars later than B3 show evidence for an abnormal 
law, in the sense that E(J - K^) is larger than expected from 
the measured E(B - V). We interpret this result as showing that 
stars later than B3 have infrared excess emission, rather than 
non-standard reddening. This excess should be a consequence 
of the presence of some remnant of the disk from which the 
stars formed and points to a very young age for NGC 1893. 

The presence of this emission excesses likely accounts for 
the existence of some dispersion in t he distance estim ates for 
NGC 18 93, ranging from DM= 13.1 dTapia et all 199 lb to 13.9 



(Paper I), as the different procedures used to deredden the data 
will treat this excess differently. In view of this difficulty, we 
prefer to assume the distance derived from spectroscopic par- 
allaxes using a standard R = 3.1 reddening law, i.e., a value 
d ~ 5 kpc, as this is in line with the dis tances to tracers of the 



Outer Arm in the vicinity of NGC 1893 (Neg ueruela & Marco 
20031) . 



Since the area covered in this in vestigati on is larger than 
the area cover ed by our photometry dPaper J) , we can use the 
photometry of iMassey et al.1 (119951) to identify B-type mem- 
bers outside the cluster core. In the B/(U - B) diagram, they 
form a very well defined sequence at bluer colours than any 
other st ar in the field. Comparison to the estimated spectral 
types of iPaper II shows that B-type members are delimited by 
(U - B) < -0.1 and B < 15.5. This agrees well with the as- 
sumption of standard reddening, as a B9 star has (U - B\) = 
-0.57 and therefore for the average reddening to the cluster 
E(U - B) = Q.12E(B - V) = 0.38, a B9 star should have 
(U - B) - -0.19. Three foreground stars can be easily iden- 
tified because of their bright B and position in the V/(B - V) di- 
agram, leaving ~ 60 stars along the cluster sequence. Of them, 
25 stars are spectroscopically confirmed to be B2 or earlier: 5 
O-type stars and 20 in the BO-2 range (of which 4 show emis- 
sion lines). 



6.2. Age 

Is there an age spread in NGC 1893? We do not find any ev- 
idence for a deviation from the ZAMS for almost any mem- 
ber. The only star that seems (slightly) evolved is the 07.5 V 
S3R1N5. However, there are other possible interpretations to 
its spectrum. For example, it could be the composite of two O- 
type stars with slightly different spectral types. We believe this 
to be rather possible, as this star lies in the densest region of the 
cluster and its spectroscopic distance modulus is rather shorter 
than the average for the cluster. Moreover, it is surrounded by 
PMS stars and shows a strong IR excess. Therefore we do not 
think that this is an evolved star, especially as nearby stars of 
earlier spectral type do not show any sign of evolution. The fact 
that an 04 V and an 05.5 V star are still close to the ZAMS 
places a strict upper limit on the a ge of the cluster at 3 My r 
and supports an even younger age (IMeynet & Maederll2003l) . 
Comparison with PMS isochrones (Fig. [6]l suggests an age of 
^ 2 Myr, in good agreement with the fact that stars as massive 
as B3-4V still show significant E(J - K&) excesses, likely 



due to the presence of remnants of the disks from which they 
formed. 

However, it is obvious that, while an important population 
of B-type stars, some of them as late as B8 - AO, is already set- 
tling into the ZAMS, some massive stars, with spectral types 
in the B1-B2 range are still in the Herbig Be phase. This 
shows that there is some spread in the formation of stars. In 
this paper, we have selected the PMS stars with the strongest 
signatures of youth (emission lines and/or strong E(H - K$) 
excess) and found that their distribution is limited to two small 
regions within the relatively large area covered by MS mem- 
bers. The fact that recent star formation is confined to these 
two regions suggests that we are observing the star formation 
process spreading from the central cluster to the neighbouring 
dark cloud. 

The lack of three-dimensional information does not allow 
us to determine the relationship of the young PMS stars close 
to the cluster centre with the members already on the main se- 
quence. Images of the area suggest that part of the dark cloud 
is partially hiding the cluster and perhaps most of the star for- 
mation is taking place on the inner wall of the cloud, which we 
cannot see. However, it seems that the population of PMS stars 
to the East of the cluster is being formed on the illuminated 
surface of the molecular cloud around the two bright pennant 
nebulae Sim 129 and Sim 130. 

6.3. Triggered star formation in Sim 130 

As described in Section [5] the area surrounding the pennant 
nebula Sim 130 contains three Herbig Be stars and several 
other emission-line stars. Other emission-line stars cover the 
area between Sim 130 and Sim 129. In the vicinit y of Sim 129, 
there are two early A-type PMS stars, S1R2N4 dPaper III) and 
SlR2N26(Table2). 

The impact of the ionising flux from the O-type stars on the 
nebulae is obvious. Their cometary aspect is due to the pres- 
ence of bright ionised fronts, taking a shape strongly resem- 
blant of a bow shock, combined with a "tail" that seems to run 
away from the centre of the cluster and is actually composed 
of bright filaments illuminated by the O-type stars. Moreover, 
similarly to other star f orming regions in the vicinity of mas- 
sive clusters (e.g., M16; lHesteretalJll996h . both nebulae show 
finger-like dust structures oriented towards the nearby O-type 
stars. 

In view of these properties, the area of recent star forma- 
tion around Sim 130 presents all the characteristics listed by 
IWalborn (2002) as typical of regions of triggered star forma- 
tion: 

- The younger (second generation) stars are associated with 
dust pillars oriented towards the O-type stars. 

- The second generation is less massive than the first. In this 
case, we have 6-7 early and mid B-type stars, as compared 
to the ~ 20 massive stars in the main cluster. 

- The more massive stars in the second generation are less 
massive than the more massive stars in the first cluster (in 
this case, the earliest spectral type around Sim 130 is Bl V, 
as opposed to 04 V in the main cluster). 
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Walbornl d2002l) finds a characteristic age difference of 
~2 Myr between the first and second generation. This may 
indeed reflect the smallest age difference that we are able to 
detect, but the observed properties of the two populations in 
NGC 1893 are not incompatible with such an age difference. 

Obviously, it may be argued that the stars around Sim 130 
may represent star formation that would have occurred regard- 
less of the effect of the nearby massive stars and that it is sim- 
ply being exposed now because of the erosion of the molecular 
cloud by the ionising flux of the O-type stars. Even though it 
is difficult to find strong arguments against this view, it should 
be noted that MSX images of the molecular cloud associated 
with IC 410 do not give any reason to suggest that star forma- 
tion is taking place anywhere except in the two areas marked 
by emission-line stars, just over the surface of the molecular 
cloud. It would be surprising if the only places in the molec- 
ular cloud in which star formation is spontaneously occurring 
happen to be, just by chance, in the process of being cleared up 
by the ionising flux of the nearby O-type stars. 

6.4. The PMS stars 

We have shown that, with the possible exception of the 
emission-line star E09, there are no obscured massive stars. We 
can consider the census of stars earlier than B3 complete. E09 
is very bright in Ks and extremely reddened. It is likely asso- 
ciated with the IRAS source 05194+33220 and therefore may 
be a massive young stellar object. 

We cannot define a strict detection limit for emission-line 
stars, but it is very unlikely that the slitless observations may 
have skipped any relatively bright (V < 15) stars. This means 
that, while there are 4 early Herbig Be stars in NGC 1893 (see 
Table 13. il l, there is one mid Herbig Be star, no late Herbig Be 
stars and only one Herbig Ae star (also, as can be seen in Fig. [5] 
the K magnitudes of all the candidate infrared excess stars are 
too faint to expect any of them to be B-type stars). This fact 
is difficult to interpret. Obviously, any reasonable IMF should 
result in a rather larger number of intermediate-mass stars than 
massive stars. Moreover, the more massive the star is, the ear- 
lier it should reach the ZAMS. Finally, the UV flux of the early 
B-type stars should help to dissolve their disks much more 
quickly than those around lower mass stars. 

How can we then explain the excess of early Herbig Be 
stars with respect to other emission-line stars? The possibility 
that early Herbig Be stars retain their disks longer than later 
type Herbig ABe stars is counterintuitive. If the star formation 
process in the two active areas is very recent, it may be pos- 
sible that only the early Herbig Be stars have emerged from 
the parental cloud and less massive stars are too faint to be 
observed in the optical. In this view, the mid and late B stars 
that we see are associated with the first generation of massive 
stars, while the second generation is now emerging from the 
parental cloud, likely because the UV photons from the O-type 
stars are photodissociating the cloud around them. The later- 



type emission line stars should then also be associated with the 
first generation of massive stars. 

6.5. The IMF 

The extreme youth of NGC 1893 offers a good prospect for 
determining the IMF of a population just emerging from the 
parental cloud. For this, deep infrared observations would be 
needed in order to probe the low-mass stellar populations. 
However, some difficulties stand out. 

First, as discussed above, it is possible that part of the clus- 
ter is obscured by parts of the dark cloud. Assuming typical 
masses for spectral types, the observed distribution of mem- 
bers is 5 stars with 25M© < M, < 60 M Q (O4-07.5), 20 
stars with 8M < M, < 16M Q (B0-B2) and ~ 35 stars with 
3M Q < M, < 8M Q (B2.5-B9). This distribution seems too 
biased towards early spectral types for a nor mal IMF. For a 
Salpeter IMF, which is valid in this mass range (lKroupal200lh . 
we would expect three times more intermediate mass (B2.5 - 
B9) than massive(<B2) stars. 

Ignoring any incompleteness due to multiplicity, the 25 
massive stars with known spectral types result in a mass ~ 
400Mq. Assuming for the sake of argument that the deficit in 
intermediate mass stars is due to observational effects and the 
IMF is standard (i.e.,|Rroupa 2001), this mass implies a total 
mass M c i ~ 22OOM for NGC 1893. This estimate is a lower 
limit, as there are reasons to believe that a substantial fraction 
of the most massive stars in the cluste r are, at least , binaries: a 
few radial velocity measurements by Ijonesl (119721) show most 
of them to display large velocity variations. 

A second factor suggesting obscuration of part of the clus- 
ter is its shape. The distribution of members is traced in Fig. [7] 
If we take the intermediate mass stars as best tracers of the clus- 
ter extent, it is difficult to assign a morphological type or even 
define a centre. The main concentration of stars appears just on 
the edge of the molecular cloud. The conspicuous absence of 
any likely member to the West of the cluster core strongly sug- 
gests that NGC 1893 is located on the back side of the molec- 
ular cloud associated with IC 410. 

An even more striking difficulty is the fact, evident in Fig. [71 
that there is a halo of high-mass stars surrounding the clus- 
ter in areas where there are essentially no intermediate-mass 
members. This is more obvious to the East of the cluster, be- 
yond Sim 130. If we consider the area lying between S2R3N35 
(RA: Q5 h 2y n \y)) and the edge of Figure UJ it contains the 8 
members identified in Fig. [T] These objects lie at distances 
of 8' to 12' from the cluster core (corresponding to 12 to 18 
pc at d ~ 5 kpc). As seen in Table 3, 7 of them have spec- 
tral types in the B0 - 1 range and one is a mi d B-type star. 



This a rea is fully covered by the photometry of iMassey et al 



3 This source is identified in SIMBAD with the 07.5 V member 
HD 242935, but its coordinates coincide much better with E09. 



dl995h . which provides only three other photometric members. 
The brightest one, fMJD951 J052325.1 3+332609.8 was clas- 
sified as B1.5 by iMassev et al. I (Tl995h . while the two other 
photometric members, [MJD95] J052336.60+332905.5 and 
J052339. 25+333839. 7 have colours and magnitudes appropri- 
ate for mid-B stars (note that J052339.25+333839.7 falls out- 
side the area covered in this investigation and lies outside 
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Fig. [7] ~ 4' to the North of the Northernmost members dis- 
played). Therefore, this area includes 8 early B stars, 3 mid 
B stars and no late B-type photometric members. How do we 
arrive at such surprising mass distribution? 

An obvious candidate for an explanatio n is dynamical elec- 
tion fr om the cluster core. As discussed by lLeonard & Duncan! 
dl990h . massive clusters containing hard binaries with two 
components of similar mass may be quite effective at eject- 
ing stars via dynamical ejections. The majority of stars ejected 
will be B-type stars and their ejection velocity will be in- 
versely proportional to their mass. These factors could explain 
a concentration of early B-type stars in the vicinity of a young 
mass ive cluster, but the effic iency at ejecting stars estimated 
bv lLeonard & Duncan! fl990) seems much lower than that re- 
quired to explain the population of objects around NGC 1893. 
Howev er, more recent work by iPflamm-Altenburg & Kroupa 
(2006) suggests that, if massive stars are mainly born as part 
of multiple s ystems, the ejection rates c an be much higher than 
estimated by lLeonard & Duncan! (11990 ). According to their re- 
sults, a cluster with a mass comparable to that of NGC 1893 
could lose ~ 75% of its high mass stars in 1 - 2 Myr. 

In this respect, it is interesting to note that there is one 
further early- t ype s tar, LS V +33°31, classified B0.5V by 



Massey et al .1 (119951) . lying about 24' away from the cluster 
core, which falls along th e sequence o f clust er members. In 
the whole area covered bv lMassev et al.l (1 1995b . there are only 
three more objects that could be photometric members, all very 
distant from the cluster (d > 15') and all compatible with be- 
ing late B-type stars [MJD95] J052425 .46+33 1544.7, [MJD95] 
J052346.82+33 1440.6 and [MJD95] J0525 16.96+332403.9. 
Even more striking is the fact that HD 242908, nominally the 
most massive star in the cluster, lies at some distance from the 
main bulk of the cluster, in an area where very few other mem- 
bers are found. Intriguingly, if we do not count stars in this 
massive star halo, the ratio between massive and intermediate- 
mass stars is close to standard. 

Radial velocity measurements of the stars in the halo of 
NGC 1893 could provide a test on this hypothesi s, as the sys- 
temic velocity of the cluster is close to zero (e.g.. |joneslll972l) 
and any measured components should be due to runaway ve- 
locities. 

One further issue to take into account is the complication 
arising from the presence of sequential star formation. If sites 
of triggered star formation (such as Sim 130) contribute only 
stars less massive than ~ 12M» (Bl V), the total integrated 
(first generation + triggered generations) population will have 
a steeper IMF than the original first generation. If we observe 
this cluster in a few Myr, there will be no way of telling which 
stars have formed at which time. Of course, this has a bearing 
on how we can define an instantaneous IMF. 



7. Conclusions 

1 . We have found a population of emission-line PMS stars in 
NGC 1893. The brightest among them cover the range from 
B 1 to late F, with an obvious overpopulation of early B-type 
stars. Emission line stars appear only in two regions of the 
cluster. 



2. We have identified a number of faint objects with high val- 
ues of (J - Ks) that seem to show an infrared excess. These 
objects concentrate around the emission-line stars, indicat- 
ing that they are also PMS stars. 

3. All the stars later than B3 show evidence for an infrared ex- 
cess, even though the main sequence is well traced down to 
AO in the optical. This infrared excess increases as we move 
to later spectral types, strongly suggesting that it arises 
from the remnant of a disk. 

4. The age of NGC 1893 is constrained to be < 3 Myr by the 
presence of main-sequence 04 and 05 stars and likely to be 
^ 2 Myr. This makes NGC 1893 one of the youngest clus- 
ters to be visible in the optical. It is very likely in the pro- 
cess of emerging from its parental cloud and perhaps more 
members lie hidden by dark portions of the cloud. If this 
is the case, they are quite faint and infrared observations 
reaching deeper than 2MASS are needed to detect them. 

5. The area around the cometary nebulae Sim 129 and 
Sim 130 shows the highest number of emission-line and 
IR-excess PMS stars. Three B1-B4 Herbig Be stars cluster 
around Sim 130. This is likely to be a region of more recent 
star formation, triggered by the ionisation front generated 
by the O-type stars. 

6. A second region containing emission-line stars and IR- 
excess PMS candidates lies on the interface between the 
cluster core and the molecular cloud. Here we could have 
another area of triggered star formation partially hidden by 
the molecular cloud. On the very edge of the cloud, we find 
the emission line object E09, which, with Ks = 9.4 and 
(J - Ks) = 3.0, could be a massive very young stellar ob- 
ject. 

7. The picture of star formation emerging from our study of 
NGC 1893 is a rather complex one, with sequential star for- 
mation resulting in several slightly non-coeval populations 
and dynamical ejection depopulating the cluster of massive 
stars at a very young age. 
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